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We present a new measurement of the unitarity triangle angle 03 using a Dalitz plot analysis 
of the KgTY + TY~ decay of the neutral D meson produced in B — > K^ 1 decays. The method 
exploits the interference between D° and D° to extract the angle 03, strong phase 8 and the ratio r 
of suppressed and allowed amplitudes. We apply this method to a 605 fb _1 data sample collected 
by the Belle experiment. The analysis uses three decays: B — > DK ± , and B ± — > D*K ± with 
D* — > Dn° and D* — > Dj, as well as the corresponding charge-conjugate modes. From a combined 
maximum likelihood fit to the three modes, we obtain 3 = 78 A°t\°il ± 3.6°(syst) ± 8.9°(model). 
CP conservation in this process is ruled out at the confidence level (1 — CL) — 5 x 1CP 4 , or 3.5 
standard deviations. 

PACS numbers: 12.15.Hh, 13.25.Hw, 14.40.Nd 



I. INTRODUCTION 

Determinations of parameters of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix [l[ are important as 
checks on the consistency of the Standard Model, and 
as ways to search for new physics. Among the angles of 
the CKM unitarity triangle, 03 (also widely known as 7) 
is the least-well constrained by direct measurements, so 
new results are of particular interest. The principal ex- 
perimental resource is CP violation in the family of de- 
cays B — > DK: various methods for extracting a 03 mea- 
surement have been proposed [l-Q , following the original 
discussion of direct CP violation measurement by Bigi, 
Carter, and Sanda The most sensitive technique re- 
lies on three-body final states 0,@1 such as Kg 7r + 7r~. 

In the Wolfenstein parameterization of the CKM ma- 
trix elements @, the weak parts of the amplitudes that 
contribute to the decay B + — > DK + are given by 
V* h V us ~ AX 3 (for the D°K+ final state) and V* b V cs ~ 
AX 3 (p + ir/) (for D°K + ); the two amplitudes interfere 
as the D° and D° mesons decay into the same final 
state Kgir + ir~. Assuming no CP asymmetry in neu- 
tral D decays, the amplitude for the process B ± — > 
(KgTr + TT~)r)K ± as a function of the Dalitz plot variables 
m\ = m 2 K o n+ and m 2 _ = m^^ is 

M± = f(m%,m^) + re^* s ^f(m%,m%), (1) 

where /(mi, m 2 _) is the amplitude of the D° — > Kgir + ir~ 
decay, r is the ratio of the magnitudes of the two inter- 
fering amplitudes, and S is the strong phase difference 



between them. The D° — > Kgir + ir~ decay amplitude / 
can be determined from a large sample of flavor-tagged 
D — > Kgir + ir~ decays produced in continuum e + e~ an- 
nihilation. Once / is known, a simultaneous fit to B + 
and B~ data allows the contributions of r, 3 and 8 to be 
separated. The method has a two- fold ambiguity: (03, S) 
and (03 + 180°, 6 + 180°) solutions cannot be separated. 
We always choose the solution with < 03 < 180°. Wc 
neglect the effects of charm mixing in this formalism. 
Given the current precision of 03 and the constraints on 
the D° mixing parameters (x n,y n ~ 0.01 [lOj]), these ef- 
fects can be safely neglected [1JJ, although it is possible 
to take them into account if they appear to be significant 
for future precision measurements. References and 
[lH give a more detailed description of the technique. 

This method can be applied to other decay modes: 
in addition to B + — > DK + , 1 excited states of neutral 
D and K can also be used, although the values of 6 
and r can differ for these decays. Both the BaBar and 
Belle collaborations have successfully applied this tech- 
nique to B^- ->• DMifMi modes with D° decaying to 
KgW + TT~ (l34l6| . In addition, the BaBar collaboration 
reported a measurement of 03 using the B^ — >• DK^ 
mode with the D decaying to the KgK + K~ [T|| and 
7r°7r + 7r~ [l7| final states. 

Here we present a measurement of 03 using the modes 



1 Charge-conjugate modes are implied throughout the paper unless 
noted otherwise. 
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B+ DK+ and B+ -> D* K+ with £) -> ^tt+tt", 
based on a 605 fb _1 data sample (657 x 10 6 Si? pairs) 
collected by the Belle detector at the KEKB asymmetric 
e + e~ collider. The Belle detector is described in detail 
elsewhere [H, HH . It is a large-solid-angle magnetic spec- 
trometer consisting of a silicon vertex detector (SVD), a 
50-layer central drift chamber (CDC) for charged particle 
tracking and specific ionization measurement (dE/dx), 
an array of aerogel threshold Chercnkov counters (ACC), 
time-of- flight scintillation counters (TOF), and an array 
of CsI(Tl) crystals for electromagnetic calorimctry (ECL) 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux return located 
outside the coil is instrumented to detect Kl mesons and 
identify muons (KLM). 

The results presented in this paper supersede our pre- 
vious measurement based on a sample of 386 x 10 6 
BB pairs [H]. In addition to B ± -)• DA' 1 and the 
B ± — > D*K ± mode with D* — > Dtt° , this analysis ex- 
ploits B ± -> D*K± with D* -> £> 7 . The D* —} 
mode has nearly the same parameters as B ± — > D*K ± 
with D* — > Dtt°, the only difference being that due to 
the opposite C parities of the 7 and ir°, the strong phases 
for these modes differ by 180° [lo||. This provides an ad- 
ditional cross-check for the analysis and allows system- 
atic uncertainties in the combined measurement to be 
reduced. The analysis procedure is also improved. It 
uses additional variables in the maximum likelihood fit 
for the separation of signal from background; this allows 
one to relax some selection requirements, thus increasing 
the sample size. 



II. EVENT SELECTION 

The decay chains B+ -> DK+ and B+ -> D*K+ with 
D* — > Dw° and D* — > Dj are selected for the analysis. 
The neutral D meson is reconstructed in the KgTr + ir~ 
final state in all cases. We also select decays D*~ —> 
D°ir~ produced via the e + e~ — ► cc continuum process as 
a high-statistics sample to determine the D° 
decay amplitude. 

Charged tracks are required to satisfy criteria based 
on the quality of the track fit and the distance from 
the interaction point. We require each track to have a 
transverse momentum greater than 100 MeV/c. (The 
reference axis is given by the direction of the e + beam.) 
Separation of kaons and pions is accomplished by com- 
bining the responses of the ACC and the TOF with the 
dE/dx measurement from the CDC. Photon candidates 
are required to have ECL energy greater than 30 MeV. 
Neutral pion candidates are formed from pairs of pho- 
tons with invariant masses in the range 120 MeV/c 2 to 
150 MeV/c 2 . Neutral kaons are reconstructed from pairs 
of oppositely charged tracks with an invariant mass M W7r 
within 7 MeV/c 2 of the nominal Kg mass, and forming 
a vertex more than 1 mm from the interaction point in 
the transverse plane. 



To determine the D° 



A°7r + 7r- 



decay amplitude 

+ - 



AgTT+TT 



we use D*^ mesons produced via the e 1 e — > cc con 
tinuum process. The flavor of the neutral D meson is 
tagged by the charge of the slow pion (which we de- 
note as tt s ) in the decay D*~ — > D°ir~. The slow pion 
track is fitted to the D° production vertex to improve 
the momentum and angular resolution of the n s . To se- 
lect neutral D candidates we require the invariant mass 
of the K s Tr + n~ system to be within 11 MeV/c 2 of the 
D° mass. To select events originating from a D*~ decay 
we impose a requirement on the difference AM of the 
invariant masses of the D*~ and the neutral D candi- 
dates: 144.9 MeV/c 2 < AM < 145.9 MeV/c 2 _ Suppres- 
sion of the combinatorial background from BB events is 
achieved by requiring the D*~ momentum in the center- 
of-mass (CM) frame to be greater than 3.0 GeV/c. The 
number of events in the signal region is 290.9 x 10 3 ; the 
background fraction is 1.0%. 

The selection of B candidates is based on the CM 
energy difference AE = ^ Ei — E\, eam and the beam- 
constrained B meson mass M bc = v^beam — Epi) 2 , 
where -©beam is the CM beam energy, and Ei and pi are 
the CM energies and momenta of the B candidate decay 
products. We select events with M bc > 5.2 GeV/c 2 and 
|Ai?| < 0.15 GeV for further analysis. We also impose a 
requirement on the invariant mass of the neutral D can- 



didate as above: \M K o 7 



M D o\ < 11 MeV/c 2 . To 



obtain the Dalitz plot variables m 2 , and m 2 _ , a kinemati- 
cal fit is employed with the constraint that the Ag7r + 7r~ 
invariant mass be equal to M D o . 

We consider two major background sources in our data: 
the continuum process e + e~ — >• qq, where the light com- 
ponent with q = u,d,s and the charmed component are 
treated separately; and BB decays, where events with 
real D° (due to B^ — > Dir^ etc.) are treated separately. 
To suppress background from continuum events, we cal- 
culate two variables that characterize the event shape. 
One is the cosine of the thrust angle cos#thr, where #thr 
is the angle between the thrust axis of the B candidate 
daughters and that of the rest of the event, calculated 
in the CM frame. The other is a Fisher discriminant T 
composed of 11 parameters 2l|: the production angle of 



the B candidate, the angle of the B thrust axis relative 
to the beam axis, and nine parameters representing the 
momentum flow in the event relative to the B thrust axis 
in the CM frame. In the first stage of the analysis, the 
(Mb c , AE) distribution is fitted in order to obtain the 
fractions of the background components, and we require 
|cos0 t hr| < 0.8 and T > -0.7. In the Dalitz plot fit, 
we do not reject events based on these variables (as in 
the previous analysis [13]), but rather use them in the 
likelihood function to better separate signal and back- 
ground events. This leads to a 7-8% improvement in the 
expected statistical error. 

The AE and M bc distributions for the B + -> DK + 
mode are shown in Fig. [T] (a), (b). For the selected 
events a two-dimensional unbinncd maximum likelihood 
fit in the variables M\, c and AE is performed, with the 
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fractions of continuum, BB and B^ 1 — > D^n^ back- 
grounds as free parameters, and their distributions fixed 
from generic MC simulation. (The continuum compo- 
nent is also split into (u, d, s) and charm components in 
the figure, based on fractions in the MC.) The resulting 
signal and background fractions are used in the Dalitz 
plot fit to obtain the event- by- event signal to background 
ratio. A more detailed description of the two-stage pro- 
cedure is given in Section IIV1 The number of events in 
the signal box (M bc > 5.27 GeV/c 2 , \AE\ < 30 MeV, 
cos0 t hr| < 0-8, J 7 > —0.7) is 756, with a signal purity 
of (70.5 ± 1.2)%. The (M bc , AE) fit yieldsa continuum 
background fraction of (17.9 ± 0.7)%, a BB background 
fraction of (7.3 ± 0.5)%, and a B ± — > Dn^ background 
fraction of (4.3 ±0.3)% in the signal box. Figure [5] shows 
the distributions of cos# t hr and J- variables in the M bc , 
AE signal region for the B + — > DK + mode. The distri- 
butions for the other modes are similar. 



To select B^ 



D*K + events with D* 



addition to the requirements described above, we require 
that the mass difference AM of neutral D* and D candi- 
dates satisfies 140 MeV/c 2 < AM < 144 MeV/c 2 . The 
AE and M bc distributions for this mode are shown in 
Fig. [T] (c), (d). The background fractions are obtained 
in the same way as for B^ — > DK mode. The number 
of events in the signal box is 149, with (79.7 ± 2.5)% 
signal purity. The continuum background fraction is 
(5.7±0.7)%, the BB background fraction is (7.6±1.9)%, 
and the B ± ->■ D*n ± background fraction is (7.0±1.3)%. 
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Selection of the B+ —> D* K+ mode with D* -> D7 
is performed in a similar way. The photon candidate is 
required to have an energy greater than 100 MeV, and 
the mass difference requirement is AM < 152 MeV/c 2 . 
Due to the larger number of background sources for this 
mode, the treatment of background differs. The BB 
background is subdivided into events with combinato- 
rial D, studied using a generic MC sample; and those 
with real neutral D mesons, for which a dedicated simu- 
lation of each component is performed. The fractions of 
background components are obtained from an unbinned 
4D fit of the distribution of variables Mb c , AE, cos^thn 
and J- . The relative fractions of BB backgrounds with 
a real D° (except for B ± ->■ D*tt ± and B ± ->■ DK±) 
are fixed according to their PDG branching ratios [22 1 
and MC efficiencies. The AE and M bc distributions for 
this mode are shown in Fig. [T] (e), (f). The number of 
events in the signal box is 141, and the signal purity 
is (41.7 ± 3.6)%. The continuum background fraction 
is (15.8 ± 1.3)%, the fraction of BB background with 
combinatorial D° is (21.3 ± 3.0)%, the contribution of 
B ± D*7T ± , IT ->■ Dj is (6.5 ± 1.2)%, and the fraction 
of the rest of BB events with real D° is (14.7 ± 1.1)%. 



The Dalitz distributions of D — > K[ 



decay in the 



FIG. 1: AE and M hc distributions for the B+ -> DK+ (a,b), 
B+ D*K+ with D* Div (c,d), and B+ -+ D*K+ with 
D* — > D7 (e,f) event samples. Points with error bars are 
the data, and the histograms are fitted contributions due to 
signal, misidentified B — > D^*'tt events, and BB, charm, 
and (u,d,s) backgrounds; in (e), a B — > DK ± contribution 
with random photon is also included. AE distributions are 
plotted with a Mb c > 5.27 MeV/c 2 requirement; M^ c distri- 
butions use a I AS I < 30 MeV requirement. |cos(?thr| < 0.8 
and T > —0.7 requirements are used in all the plots. 




signal box for each of the B ± -> DK and B ± -> D*K ± 
processes are shown in Fig. [3] 



2 4 
Fisher discr. 



FIG. 2: cos# t hr and T distributions for the B + — > DK + 
event sample. Points with error bars show the data with 
M bc > 5.27 MeV/c 2 and \AE\ < 30 MeV requirements, and 
the histograms are fitted contributions due to signal, BB, 
charm, and (u, d, s) backgrounds. 



5 




0.5 1 1.5 2 2.5 3 

m! (GeV 2 /c 4 ) 



0.5 1 1.5 2 2.5 3 

ml (GeV 2 /c 4 ) 




0.5 1 1.5 2 2.5 3 

mi (GeV 2 /c 4 ) 



0.5 1 1.5 2 2.5 3 

ml (GeV 2 /c 4 ) 



FIG. 3: Dalitz distributions of D — > K^-k + -k decays from 
selected B ± DK ± (a,b), B ± -> D*K ± with L>* ->• Dtt 
(c,d), and S — >• D* K ± with D* — >■ D7 (e,f), shown sepa- 
rately for B~ (left) and £> + (right) tags. 



III. DETERMINATION OF THE D° ArgTr+Tr" 
DECAY AMPLITUDE 

As in our previous analysis [l2[, the D° — > Kgn + TT~ 
decay amplitude is represented using the isobar model. 
The list of resonances is also the same, the only dif- 
ference being the free parameters (mass and width) of 
the K* (892) ± and p(770) states. A modified amplitude, 
where the scalar 7T7t component is described using the 
K- matrix approach [23[, is used in the estimation of the 
systematic error. 

The amplitude / for the D° — > Kgir + ir~ decay is de- 
scribed by a coherent sum of N two-body decay ampli- 
tudes and one non-resonant decay amplitude, 

N 

f(ml,rn 2 _) = ^2a j e i ^A j (ml,m 2 _) + a NR e iSNR , (2) 
3=1 

where Aj (m 2 ^ , m 2 _ ) is the matrix element, aj and £j 
are the amplitude and phase of the matrix element, 
respectively, of the j-th resonance, and qnr and £nr 



are the amplitude and phase of the non-resonant com- 
ponent. The description of the matrix elements fol- 
lows Ref. [iH • We use a set of 18 two-body am- 
These include five Cabibbo-allowed am- 
A'*(892)+tt-, A*(1410)+7r-, A *(1430)+7r-, 
and A*(1680)+7r-; their doubly Cabibbo- 



plitudes. 
plitudes: 
A|(1430)- 



suppressed partners; and eight amplitudes with Ag and 



a 7T7T resonance: 



i- s p, A 5 w, A°/ (980), A°/ 2 (1270), 
^s/o(1370), Agp(1450), K^ax and K%a 2 . 



We use an unbinned maximum likelihood technique to 
fit the Dalitz plot distribution to the model described by 
Eq.f2]with efficiency variation, background contributions 
and finite momentum resolution taken into account. The 
free parameters of the minimization are the amplitudes 
aj and phases £j of the resonances, the amplitude clnb, 
and phase ^nr of the non- resonant component, and the 
masses and widths of the a\ and a 2 scalars. We also 
allow the masses and widths of the A*(892) + and p(770) 
states to float. 

The procedures for determining the background den- 
sity, the efficiency, and the resolution are the same as in 
the previous analyses [TH, [l4[ . The background density 
for D° — > A^7r + 7r - events is extracted from AM side- 
bands. The shape of the efficiency over the Dalitz plot, 
as well as the invariant mass resolution, is extracted from 
the signal Monte-Carlo (MC) simulation. 

The fit results are given in Table fl] The fit frac- 
tion for each mode is defined as the ratio of the in- 
tegrals of the squared absolute value of the amplitude 
for that mode, and the squared absolute value of the 
total amplitude. The fit fractions do not sum up to 
unity due to interference effects. The parameters ob- 
tained for the o\ resonance (M ai = (522 ± 6) MeV/c 2 , 
r ai = (453 ± 10) MeV/c 2 ) are similar to those found 
by other experiments |25|, [26[. The second scalar term 
(T2 is introduced to account for a structure observed at 
<~ I.lGeV 2 /c 4 : the fit finds a small but signifi- 
cant contribution with M a2 = (1033 ±7) MeV/c 2 , T a2 = 
(88 ± 7) MeV/c 2 . Allowing the parameters of the dom- 
inant K* (892)+ and p(770) resonances to float results 
in a significant improvement in the fit quality. We ob- 
tain M( K* (892)) = (893.7±0.1) MeV/c 2 , T(A*(892)) = 
(48.4 ±0.2) MeV/c 2 , M{p) = (771.7 ±0.7) McV/c 2 , and 
T(p) = (136.0 ± 1.3) MeV/c 2 . 

We perform a \ 2 test using 54x54 bins in the region 
bounded by m\ = 0.3 GeV 2 /c 4 and 3.0 GeV 2 /c 4 . The 
bins with an expected population of less than 50 events 
arc combined with adjacent ones. We find y 2 /ndf = 2.35 
for 1065 degrees of freedom (ndf), which is large. We 
find that the main features of the Dalitz plot are well- 
reproduced, with some significant but numerically small 
discrepancies at peaks and dips of the distribution. In 
our final results we include a conservative contribution 
to the systematic error due to uncertainties in the D° 
decay model. 
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TABLE I: Fit results for D° 
statistical only. 
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decay. Errors are 



dimensional PDF used for the fit is thus expressed as 
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1 49 + 004 

V ' . 11./ _i_ W . W V / J 


325.4 ± 1.3 


526 ± 016 


i\~*(l410) + 7r~ 


0.65 ± 0.05 


120 ± 4 


0.49 ± 0.07 


K*(1410)-7r + 


0.42 ±0.04 


253 ±5 


0.21 ±0.03 


A' *(1430) + 7T- 


2.21 ±0.04 


358.9 ±1.1 


7.93 ±0.09 


X *(1430)-tt + 


0.36 ±0.03 


87 ±4 


0.22 ±0.04 


isT2(1430) + 7r- 


0.89 ±0.03 


314.8 ± 1.1 


1.40 ±0.06 


X 2 *(1430)-tt + 


0.23 ±0.02 


275 ±6 


0.093 ±0.014 


K*(1680) + tt- 


0.88 ±0.27 


82 ± 17 


0.06 ±0.04 


K*(1680)"7r + 


2.1 ±0.2 


130 ±6 


0.30 ±0.07 


non-resonant 


2.7 ±0.3 


160 ±5 


5.0 ± 1.0 



IV. DALITZ PLOT ANALYSIS OF B+ ->- D M K + 
DECAYS 

As in our previous analysis [ijj and in analyses car- 
ried out by the BaBar collaboration [IH [l6| , we fit the 
Dalitz distributions of the B + and B~ samples sepa- 
rately, using Cartesian parameters x± = r± cos(±03 + S) 
and y± = r± sin(±03 + S), where the indices "+" and 
"— " correspond to B + and B~ decays, respectively. In 
this approach the amplitude ratios (r+ and r_) are not 
constrained to be equal for the B + and B~ samples. Con- 
fidence intervals in r, ^3 and S are then obtained from 
the (x±,y±) using a frequentist technique. The advan- 
tage of this approach is low bias and simple distributions 
of the fitted parameters, at the price of fitting in a space 
with higher dimensionality (a;+, y+, X-, yJ) than that of 
the physical parameters {r,<f>3,$)\ see Section IVl 

Following the procedure described in Section UU back- 
ground events for the B ± -> DK ± and B ± — > D*K ± , 
D* — > Dtt q modes are classified into four components: 
e + e~ — > qq (where q = u,d,s), charm, BB (except for 
B ± _^ £)(.*) n±) an d B ± —} L>(*)tt ± background. This is 
a refinement of the previous analysis, where three back- 
ground components were used, without separation of the 
continuum background into (it, d, s) and charm. In the 
case of the B ± -> D*K ± mode with D* ->■ D7, the BB 
background is divided into events with combinatorial D, 
and seven types of events with real D mesons (including 
modes with a neutral or charged B meson decaying to 
D^*' and a K, ir or p-meson). 

The distributions of each of the background compo- 
nents are assumed to be factorized into products of a 
Dalitz plot distribution (m 2 t _,m 2 _) 1 and distributions in 
(M bc ,A£), and (cos 9^, F). The shapes of these dis- 
tributions arc extracted from MC simulation. The six- 



Pi(m 



_) Pi (M hc ,AE) Pi (cose thl ,F), (3) 



where the index i runs over all background contribu- 
tions and signal. The distributions pi(Mb c ,AE) and 
Pi(cos9thr,J~) are parameterized functions. The param- 
eterization of pi(Mb c , AE) differs for different compo- 
nents: sums of two two-dimensional Gaussian distribu- 
tions with correlations for signal and B + — > D^K + - 1 
products of the empirical shape proposed by the AR- 
GUS collaboration in M^ c 27| and a linear function in 



AE for qq, charm and BB components. In addition, the 
parameterization for BB background includes a prod- 
uct of a Gaussian peak in Mb c and a sum of exponen- 
tial and Gaussian distributions in AE. We represent 
p.;(cos#thr, J 7 ) as the product of two terms: (1) the expo- 
nential of a fourth-degree polynomial in cos# t hr, and (2) 
a sum of bifurcated Gaussian distributions in J-, where 
the mean and the various widths have a polynomial de- 
pendence on cos6>thr- The function Pi(m\,m_) is repre- 
sented by Gaussian smoothing of the MC data. 

At the first stage of the analysis (as described in Sec- 
tion |n|) we determine the relative fractions of each back- 
ground component by performing an unbinned maximum 
likelihood fit to the experimental data in M^ c and AE 
(M bc , AE, cos6» thr and 7 for B ± -> D*K ± , D* -> D7). 
The free parameters in the fit are the fractions of con- 
tinuum, BB, and B ± —> D^n^ events. The relative 
fractions of the (u, d, s) and charm components of the 
continuum background, and the relative fractions of BB 
backgrounds with real D° for B ± -> D*K ± , D* -> D7 
mode are fixed from the MC simulation. 

At the second stage, separate Dalitz distributions are 
formed for the B + and B~ samples with the signal 
requirement for M hc and AE (M hc > 5.27 GeV/c 2 , 
\AE\ < 30 McV) and no requirements for cos6*thr and 
J- . In each case, a fit with free parameters x and y is 
performed with the unbinned maximum likelihood tech- 



nique, using variables mi_, 



Mbc, AE, cos6>thr and 



J-; only the first four variables were used in the previ- 
ous analysis [l2|. Possible deviations from the factoriza- 
tion assumption for the background distribution and dis- 
agreements between MC and experimental background 
densities are treated in the systematic error. The effi- 
ciency variation as a function of the Dalitz plot variables 
is obtained from signal MC simulation and is taken into 
account in the likelihood function. 

To test the consistency of the fit, the same procedure 
as used for B + — > D^K + signal was applied to the 
B + —> D^tt + control samples. For the B ± — > 
and B ± -> D*^ (D* — > Dtt°) modes, the results are 
consistent with the expected value r ~ 0.01 for the am- 
plitude ratio. For B ± -> D*^ (D* -> D7), we find 
r = 0.05 ± 0.02, which is larger than the expected value 
by two standard deviations. Inspection of the Dalitz dis- 
tributions shows visible differences between B + and B~ 
data in this mode: we interpret the large value of r as a 
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FIG. 4: Results of signal fits with free parameters x = rcosO 
and y = rsm9 for B+ -> DK+ (a), and B+ -J- D* K+ (b) 
samples, separately for B~ and £? + data. Contours indicate 
one, two and three (for B ± DK ± ) and one standard de- 
viation regions (for _B ± — > D*K ± ) obtained in the maximum 
likelihood fit. For the B ± — > D*K ± , D* -»■ D7 mode, the sign 
of x± and y± is swapped to account for the relative strong 



phase difference of 180° with respect to the B 
D* -> Dtv° sample. 
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statistical fluctuation. 

The results of the separate B + and B~ data fits are 
shown in Fig. 2J The values of the fit parameters x± and 
y± are listed in Table [TTJ As expected, the values of x± 
and y± for the D* — > and D* — > Dn° modes from 
B ± — > D*K ± agree within the statistical errors after 
reversing the signs in one of the modes. 

V. EVALUATION OF THE STATISTICAL 
ERRORS 

We use a frequentist technique to evaluate the statis- 
tical significance of the measurements. The procedure 
is identical to that in our previous analysis [12| . This 
method requires knowledge of the probability density 
function (PDF) of the reconstructed parameters x and 
y as a function of the true parameters x and y. To ob- 
tain this PDF, we use a simplified MC simulation of the 
experiment which incorporates a maximum likelihood fit 
with the same efficiencies, resolution and backgrounds as 
used in the fit to the experimental data. 

Figure[S]shows the projections of the three-dimensional 
confidence regions onto the (r, <fo) and (<j>3,8) planes for 
the B ± ->■ DK± and B ± -> ± modes. In the results 
for the B ± ->■ D*K± mode, we combine both D* — > Dtt° 
and D* — > D-f final states, taking into account the rel- 
ative strong phase of 180° between them by swapping 
the sign of the x, y parameters for the D* — > Dj mode. 
We show the 20%, 74% and 97% confidence level regions, 
which correspond to one, two, and three standard devia- 
tions for a three-dimensional Gaussian distribution. The 
values of the parameters r, 3 and 5 obtained for the 
B ± — > DK^ and i? ± — > D* K modes separately are 
given in Table Hill The values of 4>3 in these modes agree 
within the statistical errors. In general, r and S may dif- 
fer: our results for r are similar for the two modes, while 



FIG. 5: Projections of confidence regions for the B + — > DK + 
(a,b) and B + —¥ D* K + (c,d) mode onto the (r, 03 ) and (03, ft) 
planes. Contours indicate projections of one, two and three 
standard deviation regions. 

the 8 values are distinct. 

Note that our statistical procedure gives three- 
dimensional confidence level regions. The coverage for 
the set of three parameters (r, $3,8) is exact. One- 
dimensional intervals are obtained by projecting the 
three-dimensional regions onto each of the parameter 
axes: exact coverage for this procedure is ensured only 
in the case of Gaussian errors. In our case, Gaussian 
behavior of the errors is reached when <r(r) r, and 
undercoverage (effectively, underestimation of statistical 
errors) occurs if cr(r) ~ r. The amount of undercoverage 
depends on the true value, r: errors are underestimated 
by a factor ranging from 1.4 for f = 0, to 1.03 for f equal 
to the measured value. 



VI. ESTIMATION OF SYSTEMATIC ERROR 

Experimental systematic errors come from uncertainty 
in the knowledge of the distributions used in the fit 
(i.e. Dalitz plot distributions of the background com- 
ponents, and the (Mb c , A£) and (cos#thr 5 J~) distribu- 
tions of the backgrounds and signal), fractions of differ- 
ent background components, and the distribution of the 
efficiency across the Dalitz plot. Uncertainties in back- 
ground shapes are estimated by using alternative dis- 
tributions in the fit (extracted from experimental data 
where possible). Uncertainties in the background frac- 
tions arc obtained by varying each fraction within its 
error. Possible correlations in the distributions for back- 
ground components that are not described by the formula 



8 



TABLE II: Results of the signal fits in parameters (x, y). The first error is statistical, and the second is experimental systematic 
error. Statistical correlation coefficients between x and y are also shown. Model uncertainty is not included. 



Parameter 


B+ DK+ 


B + -> D*K+, D* -> Dtt" 


B + -)■ D*K + , D* D 7 


X — 

y- 

X- — y- correlation 

x+ 

y+ 

x+ — y+ correlation 


+0.105 ±0.047 ±0.011 
±0.177 ±0.060 ±0.018 

-0.289 
-0.107 ±0.043 ±0.011 
-0.067 ±0.059 ±0.018 

±0.110 


±0.024 ±0.140 ±0.018 
-0.243 ±0.137 ±0.022 

±0.440 
±0.133 ± 0.083 ±0.018 
±0.130 ±0.120 ±0.022 

-0.101 


±0.144 ±0.208 ±0.025 
±0.196 ±0.215 ±0.037 

-0.207 
-0.006 ±0.147 ±0.025 
-0.190 ±0.177 ±0.037 

±0.080 



TABLE III: CP fit results. The first error is statistical, the second is experimental systematic, and the third is the model 
uncertainty. 



Parameter 


B + ->• DK+ mode 


B + D*K + mode 


03 

r 
8 


(80.8+^4 g ±5.0 ±8.9)° 

1R1 +0 040 , noll +0.050 

u.ioi_ 03g =c u.u±±_ 010 
(137.4+"; , ±4.0 ±22.9)° 


(73.9+ 20 2 ±4.2 ±8.9)° 
1Q6+ 00 ? 3 ± 01 3 +0 062 
(341.71 jag ± 3.2 ± 22.9)° 



(|3"| arc estimated by using independent background dis- 
tributions in the bins of Mbc A_E, cos# t i ir and T vari- 
ables. 

In case of B ± — > D*K ± , D* -> decay, an ad- 
ditional uncertainty arises from the significant cross- 
feed from the B ± D*R r± , D* Dn a mode. The 
baseline D* — > D^f fit uses x, y values obtained from 
the B ± D*#± D* L>tt° fit for modelling the 
D* — > Dtt° cross-feed; to estimate the systematic un- 
certainty, we vary x, y within their errors and also take 
x = y = 0. As an additional check, we apply a D* -± Dir° 
veto to the — > D*K ± , D* -> sample: the results 
of this fit are consistent with the baseline results within 
statistical errors. 

The procedure for estimating the uncertainty due to 
the detection efficiency is different from that in the pre- 
vious analysis [12| : here we use an alternative efficiency 
shape obtained by MC simulation from the parameter- 
ized track finding efficiency (extracted from experimental 
data) as a function of transverse momentum and polar 
angle 9. 

Compared to our previous analysis [l2j |. an additional 
source of systematic error exists due to the use of cos t hr 
and T variables in the fit. However, the use of these 
variables increases the effective signal-to-background ra- 
tio, so the total systematic error is comparable. 

Systematic errors in the physical parameters r, </>3 and 
6 are calculated from the systematic errors on the fitted 
parameters (x,y). Values (x,y) are generated according 
to Gaussian distributions with standard deviations equal 
to the corresponding total systematic errors; parameters 
r, 03 and d are then obtained for each (x,y) set, and 
the root-mean-square deviations (RMS) of the resulting 
values are calculated. We perform this procedure in two 
ways: without correlation of (x,y) biases for B + and 
B~ , and with 100% correlation between them. The larger 
RMS of the two options is chosen as the systematic error. 
The systematic errors in the x, y variables are shown in 
Table HI 



The model used for the D° — > Kgir + n~ decay ampli- 
tude is one of the main sources of error for our analysis: 
we list this contribution separately. The model uncer- 
tainty splits into two contributions: one due to imperfect 
description of the observable D° Dalitz plot distribution, 
and one due to uncertainty of the phase of the complex 
amplitude /, which is based purely on the model assump- 
tions and appears even in the case of perfect description 
of the experimental D° data. To estimate the former 
contribution, we use model variations that give a similar 
D° — > Kgir + ir~ fit quality to that of the default model. 
For the latter contribution, we take the complex phase of 
/(m^_, m?_) from models with a reduced number of reso- 
nances as in the previous analysis [l2j while keeping the 
absolute value of the amplitude the same as in the de- 
fault model. The total model uncertainty, A03 = 8.9°, 
is dominated by the uncertainty due to complex phase. 
Note that the model errors on r are highly asymmetric. 
While imperfect description of the D° density can lead 
to a bias in both directions, a wrong complex phase in- 
troduces a bias only to lower values. 

Our estimate of the model uncertainty can be consid- 
ered conservative. When the various 5*- wave terms - 
the most theoretically controversial part of the model — 
are replaced by a i^-matrix amplitude (23j . the change 
in 03 from the baseline fit does not exceed 2°. However, 
we retain our default 8.9° uncertainty as the if-matrix 
describes only part of the amplitude. 

Using a different approach, it is possible to remove 
the current model uncertainty, exploiting constraints on 
the complex phase in the D° — > Kgir + ir~ amplitude 
that can be obtained experimentally from the analysis of 
■0(3770) — > D°D° decays. Such a measurement was re- 
cently performed by CLEO [28|- The results show good 
agreement with the isobar model, however a quantitative 
estimate of the model uncertainty for a mo del- dependent 
fit is hard to obtain from these data. Instead, a model- 
independent analysis @, HI; H3] involving a binned fit of 
the D° — s- Kgir + ir~ Dalitz distribution is possible. The 
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model error in this analysis will be replaced by a statisti- 
cal error of about 1-2° due to the finite 0(3770) -> D°D° 
sample, while the statistical error associated with the B 
data sample should increase by 10-20% due to the binned 
fit procedure. At the current level of precision, this will 
not result in a significant improvement in the precision of 
03, but future analyses with larger samples of B decays 
should benefit from the model-independent technique. 

VII. COMBINED 3 MEASUREMENT 

The two event samples, B + — > DK + and B + — > 
D*K + , are combined in order to improve the sensitiv- 
ity to 03. The confidence levels for the combination of 
the two modes are obtained using the same frequcntist 
technique as for the single mode, with the PDF of the 
two measurements being the product of the probability 
densities for the individual modes. Confidence intervals 
for the combined measurement together with systematic 
and model errors are shown in Table IIVI The statistical 
confidence level of CP violation is 1 — CL = 1.5 x 10~ 4 , or 
3.8 standard deviations. With the systematic and model 
errors taken into account, CP conservation is ruled out at 
the confidence level 5 x 10 -4 , or 3.5 standard deviations. 
The systematic errors are assumed to be uncorrelated in 
this calculation; the resulting estimate is conservative, as 
most of the systematic biases are correlated between B + 
and B~ samples and thus do not introduce CP violation. 

VIII. CONCLUSION 

We report the results of a measurement of the unitarity 
triangle angle 03, using a method based on Dalitz plot 
analysis of TP — > Kgir + Tr^ decay in the process B + — > 
K + . A new measurement of 03 using this technique 
was performed based on 605 fb _1 of data collected by the 
Belle detector: 70% larger than the previous sample [l2j]. 
The statistical sensitivity of the measurement has also 
been improved by modifications to the event selection 
and fit procedure, and by adding the sample with D* 
decaying to the Dj final state. 

From the combination of B + — > DK + and 
B + — > D*K + modes, we obtain the value 03 = 
78.4°l^Jo(stat)±3.6 (syst)±8.9°(model); of two possi- 
ble solutions we have chosen the one with < 03 < 180°. 
We also obtain values of the amplitude ratios run = 
0.160lg;^(stat)±0.011(syst)+o ^o(model), and r D , K = 

0.196t^g(stat) ± 0.012(syst)+°;^2 (model). The CP 



conservation in the combined measurement is ruled out 
at the confidence level 5 x 10 -4 , or 3.5 standard devia- 
tions. 

The statistical precision of the 03 measurement is al- 
ready comparable to the estimated model uncertainty. 
However, it is possible to eliminate this model uncer- 
tainty using constraints on the D° — > Kgir + ir~ decay 
amplitude obtained by the CLEO collaboration in the 
analysis of 0(3770) D°D° decays (US]. The sta- 
tistical errors in the proposed binned fit procedure are 
10-20% larger, but the model uncertainty is replaced 
by a small (1 — 2°) statistical error due to the finite 
-0(3770) — > D°D° sample. This should result in an im- 
provement of the 03 precision in future high-statistics 
analyses. 



Acknowledgments 

We thank the KEKB group for the excellent opera- 
tion of the accelerator, the KEK cryogenics group for the 
efficient operation of the solenoid, and the KEK com- 
puter group and the National Institute of Informatics 
for valuable computing and SINET3 network support. 
We acknowledge support from the Ministry of Educa- 
tion, Culture, Sports, Science, and Technology (MEXT) 
of Japan, the Japan Society for the Promotion of Sci- 
ence (JSPS), and the Tau-Lepton Physics Research Cen- 
ter of Nagoya University; the Australian Research Coun- 
cil and the Australian Department of Industry, Innova- 
tion, Science and Research; the National Natural Sci- 
ence Foundation of China under contract No. 10575109, 
10775142, 10875115 and 10825524; the Ministry of Ed- 
ucation, Youth and Sports of the Czech Republic un- 
der contract No. LA10033; the Department of Science 
and Technology of India; the BK21 and WCU program 
of the Ministry Education Science and Technology, Na- 
tional Research Foundation of Korea, and NSDC of the 
Korea Institute of Science and Technology Information; 
the Polish Ministry of Science and Higher Education; the 
Ministry of Education and Science of the Russian Federa- 
tion and the Russian Federal Agency for Atomic Energy; 
the Slovenian Research Agency; the Swiss National Sci- 
ence Foundation; the National Science Council and the 
Ministry of Education of Taiwan; and the U.S. Depart- 
ment of Energy. This work is supported by a Grant- 
in-Aid from MEXT for Science Research in a Priority 
Area ("New Development of Flavor Physics"), and from 
JSPS for Creative Scientific Research ( "Evolution of Tau- 
lepton Physics"). 



[1] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, (1991). 

652 (1973); N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963). [3] I. Dunietz, Phys. Lett. B270, 75 (1991). 

[2] M. Gronau and D. London, Phys. Lett. B253, 483 [4] D. Atwood, G. Eilam, M. Gronau and A. Soni, Phys. 

(1991); M. Gronau and D. Wyler, Phys. Lett. B265, 172 Lett. B341, 372 (1995). 



10 



TABLE IV: Results of the combination of the B + -> DK + and B + -»■ D*K + modes. 



Parameter 


la interval 


2cr interval 


Systematic error 


Model uncertainty 


03 


(78.41%*)° 

o.i6ot°;°« 

0.196±g;g« 


54.2° < 03 < 100.5° 
0.084 < r DK < 0.239 
0.061 < r D * K < 0.271 


3.6° 
0.011 
0.012 


8.9° 

+0.050 
-0.010 
+0.062 
-0.012 


Sdk 


(136.7t£°)° 


102.2° < S DK < 162.3° 


4.0° 


22.9° 


Sd-k 


(341.9111:2)° 


296.5° < 5d*k < 382.7° 


3.0° 


22.9° 



[5] D. Atwood, I. Dunietz and A. Soni, Phys. Rev. Lett. 78, [18 

3257 (1997); D. Atwood, I. Dunietz and A. Soni, Phys. 

Rev. D 63, 036005 (2001). [19 
[6] I. I. Bigi and A. I. Sanda, Phys. Lett. B211, 213 (1988); 

A. B. Carter and A. I. Sanda, Phys. Rev. Lett 45, 952 [20 

(1980). 

[7] A. Giri, Yu. Grossman, A. Soffer, J. Zupan, Phys. Rev. [21 

D 68, 054018 (2003). 
[8] A. Bondar. Proceedings of BINP Special Analysis Meet- [22 

ing on Dalitz Analysis, 24-26 Sep. 2002, unpublished. 
[9] L. Wolfenstein, Phys. Rev. Lett. 51, 19 45 (1983). [23 
[10] A. J. Schwartz (for HFAG charm group). larXiv:0911. 14641 

[hep-ex]. [24 
[11] Y. Grossman, A. Soffer, J. Zupan, Phys. Rev. D 72, 

031501 (2005). [25 
[12] Belle Collaboration, A. Poluektov et al., Phys. Rev. D 

73, 112009 (2006). 
[13] Belle Collaboration , K. Abe et al, [26 

|arXiv:hep-ex/0308043] 
[14] Belle Collaboration, A. Poluektov et al. Phys. Rev. D [27 

70, 072003 (2004). 
[15] BaBar Collaboration, B. Aubert et al., Phys. Rev. Lett. [28 

95, 121802 (2005). 
[16] BaBar Collaboration, B. Aubert et al., Phys. Rev. D 78, [29 

034023 (2008). [30 
[17] BaBar Collaboration, B. Aubert et al., Phys. Rev. Lett. 

99, 251801 (2007). 



Belle Collaboration, A. Abashian et al., Nucl. Instr. and 
Meth. A 479, 117 (2002). 

Y. Ushiroda (Belle SVD2 Group), Nucl. Instr. and Meth. 
A 511, 6 (2003). 

A. Bondar and T. Gershon, Phys. Rev. D 70, 091503 
(2004). 

CLEO Collaboration, D. M. Asner et al, Phys. Rev. D 
53, 1039 (1996). 

Particle Data Group, C. Amsler et al, Phys. Lett. B 667, 
1 (2008). 

V. V. Anisovich and A. V. Sarantsev, Eur. Phys. J. A 
16, 229 (2003). 

CLEO Collaboration, S. Kopp et al, Phys. Rev. D 63, 
092001 (2001). 

CLEO Collaboration, H. Muramatsu et al., Phys. Rev. 
Lett. 89, 251802 (2002), Erratum- ibid: 90, 059901 
(2003). 

E791 Collaboration, E. M. Aitala et al, Phys. Rev. Lett. 
86, 770 (2001). 

ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 
241, 278 (1990). 

CLEO Collaboration, R.A. Briere et al., Phys. Rev. D 
80, 032002 (2009). 

A. Bondar, A. Poluektov, Eur. Phys. J. C 47, 347 (2006). 
A. Bondar, A. Poluektov, Eur. Phys. J. C 55, 51 (2008). 



